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Abstract 

A fundamental result relevant to spin chains and two-dimensional disordered systems is 
that the sphere sigma model with instanton coupling 9 = tt has a non-trivial low-energy 
fixed point and a gapless spectrum. This result is extended to two series of sigma models 
with 9 = tt: the SU(N)/SO{N) sigma models flow to the SU(N)x WZW theory, while the 
0(2N)/0(N) x O(N) models flow to 0(2N) X (2N free Majorana fermions). These models 
are integrable, and the exact quasiparticle spectra and S matrices are found. One interesting 
feature is that charges fractionalize when 9 = tt. I compute the energy in a background field, 
and verify that the perturbative expansions for 9 = and tt are the same as they must be. 
I discuss the flows between the two sequences of models, and also argue that the analogous 
sigma models with Sp(2N) symmetry, the Sp(2N)/U(N) models, flow to Sp(2N)i. 

1 Introduction 

One very interesting property of field theories is that they can have critical points which are 
completely unseen in standard weak-coupling perturbation theory. To prove the existence of 
these hxed points requires using alternative perturbative methods (such as the large- A r or large- 
spin expansion), or relying on non-perturbative methods. One of the remarkable features of two- 
dimensional held theories is that in some cases these non-trivial critical points can be understood 
non-pert urbatively. 

A famous example arose in the study of integer and half-integer spin chains Q. The spin- 1/2 
Heisenberg quantum spin chain is exactly solvable, and from Bethe's exact solution it is known 
that the spectrum is gapless. The obvious guess for the held theory describing the spin chain in 
the continuum limit is the sphere sigma model. This sigma model is an ,S'[/(2)-symmetric held 
theory where the held takes values on a two-sphere. The two-sphere can be parametrized by 
three fields (wi,^,^) obeying the constraint {v\) 2 + (t>2) 2 + (V3) 2 = 1. The Euclidean action is 

3 

-J2[ dxdy {d^Vif (1) 

s 1=1 ^ 

This held theory, however, is not the correct continuum description of the spin chain, because 
it is also exactly solvable, and its spectrum is gapped [g]. It was proposed that that when 
this sigma model is modified by adding an extra term called the theta term, the spectrum can 
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become gapless [jjj. The theta term is inherently non-perturbative: it does not change the beta 
function derived near the trivial fixed point at all. Nevertheless, the physics can be dramatically 
different from when 9 = ||]. In two dimensions, this can result in a non-trivial fixed point. Not 
only does the spectrum become gapless at the non-trivial critical point, but in addition, charge 
fractionalization occurs: the charges of the quasiparticles of the theory are fractions of the charge 
of the fields in the action. 

To define the theta term, field configurations are required to go to a constant at spatial 
infinity, so the spatial coordinates (x, y) are effectively that of a sphere. Since the field takes 
values on a sphere as well, the field is therefore a map from the sphere to a sphere. An important 
characteristic of such maps is that they can have non-trivial topology: they cannot necessarily 
be continuously deformed to the identity map. This is analogous to what happens when a circle 
is mapped to a circle (i.e. a rubber band wrapped around a pole): you can do this an integer 
number of times called the winding number (a negative winding number corresponds to flipping 
the rubber band upside down). It is the same thing for a sphere: a sphere can be wrapped around 
a sphere an integer number of times. An example of winding number 1 is the isomorphism from a 
point on the spatial sphere to the same point on the field sphere. The identity map has winding 
number 0: it is the map from every point on the spatial sphere to a single point on the field 
sphere, e.g. (vi(x, y),V2(x, y),v$(x, y)) = (1,0,0). Field configurations with non-zero winding 
number are usually called instantons. The name comes from viewing one of the directions as 
time (in our case, one would think of say x as space and y as Euclidean time). Since instanton 
configurations fall off to a constant at y = ±oo, the instanton describes a process local in time and 
hence "instant". Therefore, the field configurations in the sphere sigma model can be classified 
by an integer n. This allows a term 

Sg = ind 

to be added to the action, where 6 is an arbitrary parameter. Since n is an integer, the physics is 
periodic under shifts of 2tt in 9. Haldane argued that when 9 = ir, the sphere sigma model flows 
to a non-trivial critical point ||]. 

Around the same time, a similar proposal arose in some more general sigma models Q. These 
models are most easily formulated by having the field take values in the coset space G/H, where 
G and H are Lie groups, with H a subgroup of G. In this language, the two-sphere is equivalent 
to 0(3)/0(2): while the vector i>2, ^3) can be rotated by the 0(3) symmetry group, it is 
invariant under the 0{2) subgroup consisting of rotations around its axis. Thus the space of 
distinct three-dimensional fixed- length vectors (the sphere) is the coset 0(3)/0(2). In order to 
describe two-dimensional non-interacting electrons with disorder and a strong transverse magnetic 
field, one takes G = U(2N) and H = U(N) x U(N) For N = 1, one recovers the sphere sigma 
model. Pruisken conjectured that in the replica limit N — > 0, this sigma model has a critical 
point at 9 = tt. This critical point possibly describes the transitions between integer quantum 
Hall plateaus. Subsequently, a number of conjectures for non-trivial fixed points in sigma models 
have been made. For a survey of the applications of these conjectures to disordered systems, see 

The purpose of this paper is to generalize this result and prove that several infinite hierarchies 
of sigma models have non-trivial fixed points when 9 = tt. A useful tool is to study the spectrum 
and scattering matrix of the particles in an equivalent 1 + 1-dimensional formulation. I will show 
that when these sigma models have 9 = 0, they are gapped, with the spectrum consisting of 
massive particles in the symmetric representation of SU(N) and 0(2P) respectively (plus bound 
states in more general representations, and kinks with fractional charge in the latter case)). 
When 9 = tt, the spectrum consists of gapless quasiparticles which are in the fundamental 



2 



representations (vector, antisymmetric tensor,. . . ) of SU(N) and 0{2P) (the latter including 
kinks in the spinor representation). 

The stable low-energy fixed points for the sphere sigma model and all the sigma models 
discussed in this paper are Wess-Zumino-Witten (WZW) models. The WZW model for a group 
H is a sigma model where the field takes values in H x H/H ~ H. It also has an extra term, 
called the Wess-Zumino term, which has an integer coefficient k. In two dimensions, when k ^ 0, 
the model has a stable low-energy fixed point ||. The conformal field theory describing this fixed 
point is called the WZW model. 

One argument for the flow in the sphere sigma model at 9 = tt goes as follows @ . First one uses 
Zamolodchikov's c-theorem, which makes precise the notion that as one follows renormalization 
group flows, the number of degrees of freedom goes down. Zamolodchikov shows that there 
is a quantity c associated with any two-dimensional unitary field theory such that c must not 
increase along a flow. At a critical point, c is the central charge of the corresponding conformal 
field theory ||. At the trivial fixed point of a sigma model where the manifold is flat, the central 
charge is the number of coordinates of the manifold. For the sphere, this means that c = 2 at the 
trivial fixed point. This if the sphere sigma model flows to a non-trivial fixed point for 9 = tt, 
this fixed point must have 0(3) ~ SU(2) symmetry and must have central charge less than 2. 
The only such unitary conformal field theories are SU(2)k for k < 4. (The central charge of 
SU{2)^ is 3k/(k + 2); in general, the central charge of iT& is k dim(H)/(k + h), where h is the 
dual Coxeter number of H). One can use the techniques of Q to show that there are relevant 
operators at these fixed points, and at k = 2 or 3, no symmetry of the sphere sigma model 
prevents these relevant operators from being added to the action |7|]. So while it is conceivable 
that the sphere sigma model with 9 = tt could flow near to these fixed points, these relevant 
operators would presumably appear in the action and cause a flow away. However, there is only 
one relevant operator (or more precisely, a multiplet corresponding to the WZW field w itself) 
for the SU(2)i theory. The sigma model has a discrete symmetry (vx,V2,Vs) — > (— vi, — 1>2, — V3) 
when 9 = or tt; the winding number n goes to — n under this symmetry, but 9 = tt and 9 = —tt 
are equivalent because of the periodicity in 9. This discrete symmetry of the sigma model turns 
into the symmetry w — > —w of the WZW model. While the operator tr w is SU(2) invariant, it is 
not invariant under this discrete symmetry. Therefore, this operator is forbidden from appearing 
in the effective action. The operator (tru;) 2 is irrelevant, so it is consistent for the sphere sigma 
model at 9 = tt to have the SU(2)i WZW model as its low-energy fixed point. A variety of 
arguments involving the spin chain strongly support this conjecture [Q, [jjj. 

An important question is therefore whether the existence of these non-perturbative fixed 
points in sigma models at 9 = tt can be established definitively. The fact that the sphere sigma 
model has a non-trivial fixed point at 9 = it was proven in |n], 0|. This proof does not involve 
the spin-1/2 chain which motivated the result. Rather, it is a statement about the non-trivial 
fixed point in the sphere sigma model at 9 = tt. This proof utilizes the integrability of the sphere 
sigma model at 9 = and tt. Integrability means that there are an infinite number of conserved 
currents which allow one to find exactly the spectrum of quasiparticles and their scattering matrix 
in the corresponding 1 + 1 dimensional field theory. The quasiparticles for 9 = are gapped and 
form a triplet under the SU{2) symmetry, while for 9 = tt they are gapless, and form SU (2) 
doublets (left- and right-moving) [||, |i~0[ . This is a beautiful example of charge fractionalization: 
the fields (v 1, t>2, V3) form a triplet under the SU(2) symmetry, but when 9 = tt the excitations of 
the system are doublets. To prove that this is the correct particle spectrum, first one computes 
a scattering matrix for these particles which is consistent with all the symmetries of the theory. 
From the exact S matrix, the c function can be computed. It was found that at high energy c 
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indeed is 2 as it should be at the trivial fixed point, while c = 1 as it should be at the SU(2)i 
low-energy fixed point [10]. As an even more detailed check, the free energy at zero temperature 
in the presence of a magnetic field was computed for both 6 = Il2] and tt pl| . The results can 
be expanded in a series around the trivial fixed point. One can identify the ordinary perturbative 
contributions to this series, and finds that they are the same for 6 = and tt, even though the 
particles and S matrices are completely different [O]. This is as it must be: instantons and the 
6 term are a boundary effect and hence cannot be seen in ordinary perturbation theory. One 
can also identify the non-perturbative contributions to these series, and see that they differ. Far 
away from the trivial fixed point, non-perturbative contributions can dominate which allow a 
non-trivial fixed point to appear when 6 = tt even though there is none at 6 = 0. 

The purpose of this paper is to generalize these computations to the SU(N)/SO(N) and 
0(2P)/0(P) x O(P) sigma models. To have any hope of being able to take the replica limit 
N — > 0, one needs a solution for any N. The SU(N)/SO(N) sigma model reduces to the sphere 
sigma model when N = 2, while 0(2P)/0(P) x O(P) reduces to two copies of the sphere when 
P = 2. Thus these sigma models are the generalizations of the sphere sigma model with SU(N) 
and 0(2P) symmetry respectively. One difference, however, is that in general they do not allow 
a continuous 6 parameter, but instead 6 can only be zero or tt. I will show that when 6 = tt, the 
SU(N)/SO(N) and 0(2P)/0(P) x O(P) sigma models have stable low-energy fixed points. The 
corresponding conformal field theories are the SU(N)i and 0{2P)\ WZW models, respectively. 

In section 2 I define the sigma models. In section 3, I review S matrices in an integrable 
model with a global symmetry G. In section 4, I discuss the Gross-Neveu model. This is not 
a sigma model, but is closely related. These results are crucial in what follows. In section 5, 
I will find the particles and their scattering matrices for the sigma models. In section 6, I will 
do a substantial check on this picture by using the exact S matrix to compute the energy in 
the presence of a background field. In particular, I check that the perturbative contributions to 
this energy are in agreement with conventional perturbation theory, and that the expansions for 
6 = and tt are the same (but the non-perturbative contributions differ). In section 7, I discuss 
the low-energy fixed points at 6 = tt, and flows between them. I also present a conjecture for 
the low-energy behavior of the Sp(2N)/U(N) sigma model. This last section is written to be 
reasonably self-contained, so readers interested in the results but not the details (and willing to 
trust me) can skip sections 3-6. 



2 The models 

The sigma models discussed in this paper all can be written conveniently in terms of an matrix 
field <I> with action 

S = -tr J d 2 x W&dpQ (2) 

along with the constraints 



det(§) 



I 

±1 



(3) 



where / is the identity matrix. In other words, $ is always a unitary matrix with determinant 
±1, with possibly additional constraints. The constraints (|[) can easily be obtained from theories 
without constraints by adding a potentials like Atr (<^<3? — I) 2 . When A gets large, one recovers 
the restrictions (0). In theories with interacting fermions, this often results from introducing a 
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bosonic field to replace four-fermion interaction terms with Yukawa terms (interactions between 
a boson and two fermions). Integrating out the fermions then results in such potentials and hence 
the sigma model. 

The sigma models discussed in this paper are obtained by putting additional restrictions of 
the matrix field <£. The two series of models studied both require that $ be symmetric as well as 
unitary, so the matrices do not form a group: multiplying symmetric matrices does not necessarily 
give a symmetric matrix. Instead, these spaces is of G/H form. The SU(N)/SO(N) sigma model 
is obtained by requiring that $ be a symmetric N x N unitary matrix with determinant 1. The 
0(2P) /O(P) x O(P) sigma model is obtained by requiring that $ be a real, symmetric, orthogonal 
and traceless 2P x 2P matrix. In general, in two dimensions a G/H sigma model has a global 
symmetry G. The action @) and restriction (||) are invariant under the symmetry 

$ -» U$U T , (4) 

where U is a unitary matrix with determinant one. This is the most general symmetry which 
keeps $ symmetric and unitary. In the 0{2P)/0{P) x O(P) sigma models, U must be real as 
well, so G = 0(2P). The field $ in this case can diagonalized with an orthogonal matrix U, so 

$ = UAU T $ £ 0(2P)/0(P) x O(P), 

where U is in 0(2P), and A is the matrix with P values +1 and P values —1 on the diagonal. 
Different U can result in the same the subgroup leaving $ invariant is H = O(P) x O(P). 
Similarly, field configurations in the SU(N)/SO(N) sigma model can be written in the form 

<$> = UU T $ G SU(N)/SO(N) 

where U is in SU(N). The subgroup H leaving $ invariant is SO(N). For example, $ = I for 
any real U in SU(N), i.e. if U is in the real subgroup SO(N) of SU(N). This is why H = SO{N) 
here. 

These spaces G/H studied here are examples of symmetric spaces. A symmetric space G/H 
has H a maximal subgroup of G (no subgroup other than G itself contains H). The important 
property of a sigma model on a symmetric space is that it contains only one coupling constant g 
in the action. In other words, the space G/H preserves its "shape" under renormalization, with 
only the overall volume changing. The one and two-loop beta functions for all two-dimensional 
sigma models are universal and can be expressed in terms of the curvature of the field manifold 



[13]. The effect of renormalization is to increase the curvature (increase g). Even though naively 
there is no mass scale in the theory (g is dimensionless) , physical quantities depend on a scale (e.g. 
a lattice length or a momentum cutoff) as a result of short-distance effects. The renormalized 
coupling depends on this scale. At short distances, g is small, so the theory is effectively free: 
it is a theory of dimG— dimH free fields. At longer distances the theory is interacting: the 
dimG— dinxff fields interact because of the constraint (||). In renormalization-group language, 
there is an unstable trivial fixed point at g = 0, and no non-trivial low-energy fixed point unless 
one includes a theta or WZW term. 

Both sets of sigma models allow a theta term to be added to the action. This is a result long 
ago proven by mathematicians (for a discussion accessible to physicists, see ]14|]). In mathematical 
language, the question is whether the second homotopy group ^(G/H) is non-trivial. The second 
homotopy group is just the group of winding numbers of maps from the sphere to G/H, so for 
the sphere it is the integers. The general answer is that ^(G/H) is the kernel of the embedding 
of tti(H) into 7Ti(G), where %\{H) is the group of winding numbers for maps of the circle into 
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H. The rubber band on a pole example means that iti{H) is the integers when H is the circle 
= U{1) = SO(2). The only simple Lie group H for which ix\ is nonzero is SO(N), where 
tti(SO(N)) = Z2 for N > 3 and Z for N=2. Thus there are models with integer winding 
number, some with just winding number or 1, and some with no instantons at all. Integer 
winding number means that 9 is continuous and periodic, while a winding number of or 1 means 
that 9 is just or tt (just think of 9 as being the Fourier partner of n). The SU(N)/ 'SO(N) and 
0{2P)/0{P) x O(P) sigma models therefore have instantons with Z2 winding number, and so 9 
can be zero or tt. 

3 Generalities on exact S matrices 

Even though the sigma models are originally defined in two-dimensional Euclidean space (x,y), 
it is very convenient to continue to real time t = iy. The reason is to treat the field theory as 
a 1 + 1 dimensional particle theory. A fundamental property of many field theories is that all 
states of the theory can be written in terms of particles. In other words, the space of states 
is combination of one-particle states, usually called a Fock space. If the theory has a stable 
non-trivial fixed point at low energy, then the particles should be massless: the energy is linearly 
related to the momentum: E = \P\. If there is no fixed point, then the particles are massive: 
E = yJP 2 + M 2 in a Lorentz-invariant theory. In condensed-matter physics, the particles are 
often called "quasiparticles" , to emphasize the fact that the particles may not be the same as 
the underlying degrees of freedom: just because a system is made up of electrons does not mean 
that the collective excitations are electrons, or even resemble them. This phenomenon will occur 
frequently in this paper. 

One finds the scattering matrix of an integrable model by utilizing a variety of constraints: 
unitarity, crossing symmetry, global symmetries, the consistency of bound states, and the fac- 
torization (Yang-Baxter) equations. It is convenient to write the momentum and energy of a 
particle in terms of its rapidity 9, defined by E = m cosh 9, P = msinh#. Lorentz invariance 
requires that the two-particle S matrix depends only on the rapidity difference 9\ — #2 of the two 
particles. 

The invariance of the G/H sigma model under the Lie-group symmetry G requires that the 
S matrices commute with all group elements. The S matrix can then be conveniently written 
in terms of projection operators. A projection operator Vk maps the tensor product of two 
representations onto an irreducible representation labelled by k. By definition, these operators 
satisfy VkPi = SkfPk- Requiring invariance under G means that the S matrix for a particle in 
the representation a with one in a representation b means that the S matrix is of the form 



where 9 = 9 a — 9^ is the difference of the rapidities, and the are as of yet unknown functions. 
The sum on the right-hand side is over all representations k which appear in the tensor product 
of a and b; of course ^2uVk = 1- hi an integrable theory, the functions f% b (9) are determined 
up to an overall function by requiring that they satisfy the Yang-Baxter equation. This stems 
from the requirement that the S'-matrix be factorizable: the multiparticle scattering amplitudes 
factorize into a product of two-particle ones. There are two possible ways of factorizing the 
three-particle amplitude into two-particle ones; the requirement that they give the same answer 
is the Yang-Baxter equation. There have been hundreds of papers discussing how to solve this 




(5) 



k 
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equation, so I will not review this here. For a detailed discussion relevant to the sigma models 
here, see e.g. [15, 16, [17]]. Solutions arising in the sigma models will be given below. 



To obtain the overall function not given by the Yang-Baxter equation, one needs to require 
that the S matrix be unitary, and that it obey crossing symmetry. With the standard assump- 
tion that the amplitude is real for 9 imaginary, the unitarity relation S'(9)S(9) = I implies 
S(9)S(—9) = I. The latter is more useful because it is a functional relation which can be con- 
tinued throughout the complex 9 plane. Crossing symmetry is familiar from field theory, where 
rotating Feynman diagrams by 90° relates scattering of particles a, and bj to the scattering of 
the antiparticle o« with bj. In S matrix form, it says that 

S ab (in-9) = C a S ab (9)C a , (6) 

where C a is the charge-conjugation operator acting on the states in representation a. 

Multiplying any S matrix by function F(9) which satisfies F(9)F(—9) = 1 and F(m — 9) = 
F(9) will give an S matrix still obeying the Yang-Baxter equation, crossing and unitarity (this 
is called the CDD ambiguity). To determine F{9) uniquely, one ultimately needs to verify that 
the S matrix gives the correct results for the free energy. How to do this will be described in 
the following sections. Before doing this calculation, one must make sure an additional criterion 
holds: the poles of the S matrix are consistent with the bound-state spectrum of the theory. The 
S matrices of bound states are related to those of the constituents by the bootstrap relation, 
which can be formulated as follows fl5| , |i~6| . Poles of S ab matrix elements at some value 9 = 9 ab 
with 9 a b imaginary and in the "physical strip" < Im(9 ab ) < ir are usually associated with 
bound states. Each of the functions f^ b {9) has a residue at this pole Rk- Then there are bound 
states in representation k if R^k < 0, where rjk = ±1 is the parity of states in representation k 



18] (poles which do not correspond to bound states give bound states in the process obtained 





by crossing). These bound states (ab) have mass 

(^(ab)) 2 = (w-a) 2 + (m b ) 2 + 2m a m b cosh 9 ab 

The bound states are not in a irreducible representation of G if more than one of the residues 
Rk is non-zero. Then the S matrix for scattering the bound state (ab) from another particle c is 
given by 

S {ab)c = ( J] s/\Rk~\r k ) S ac (9 + i^)S bc (9 -i^)[Y: -^Vk ) , (7) 



where the projection operators Vk act on the states in representations a, b. Note that the matrices 
in (0) are not all acting on the same spaces, so this relation is to be understood as multiplying 
the appropriate elements (not matrix multiplication). 

All of the above considerations apply to massless particles as well, with a few modifications 
and generalizations [|irj, |l9| . The theories being studied are along a flow into their stable low- 
energy fixed point. Even though the quasiparticles are gapless, there is still a mass scale M 
describing the crossover: M — > is the unstable high-energy fixed point, while M — > oo is the 
stable low-energy one. The rapidity variable for a massless particle is then defined via E = me 9 , 
P = me 9 for right movers and E = me~ e , P = —me~ 9 for left movers, where m is not really 
the mass of the particle but is proportional to M (different particles can and do have different 
values of m). The S matrix still depends only on rapidity differences. Note in particular that 
the S matrices for two left movers Sll or two right movers Srr depend on the ratio E\jEi of 
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the two particles, and so do not depend on the mass scale M at all. These are determined solely 
by properties of the low-energy fixed point. On the other hand, Slr depends on the physics of 
how one flows into this fixed point. Massless particles still can have bound states like massive 
ones. These bound states show up as poles in Sll and Srr, but it does not seem possible to 
have poles in Slr in the physical strip |lC[| . 



4 The Gross-Neveu models 

Before discussing the particle content of the sigma models, it is useful to first discuss the Gross- 
Neveu models. These are not sigma models, but are closely related: they are asymptotically-free 
field theories with a global symmetry G. For groups O(N) and U(N) they were originally 
formulated in terms of fermions with a four-fermion interaction. For O(N), the action was first 
written in terms of N Majorana (real) fermions with a four-fermion interaction: 2C] 

N ( N \ 2 

S = J E faidLii + ^r^r) + 9 I E Vi 4t I ■ (8) 

For SU(N) the action can also be written in terms of fermions (this is sometimes known as the 
chiral Gross-Neveu model |^2j ) , but for general G the Gross-Neveu models are most easily defined 
in terms of the G Wess-Zumino-Witten model at level 1. The WZW model has a local symmetry 
Gl x Gr generated by chiral currents j£ and j R (the explicit definitions are given in section 7 
below). The action of the G Gross-Neveu model is then a perturbation of the G\ WZW action: 

S = S(G 1 )+gJ2 [ fLJR (9) 

The coupling g is naively dimensionless, but it has lowest-order beta function 

13(g) oc g 2 , 

where the constant of proportionality is negative. This means that for positive g the coupling is 
marginally relevant |2(|. In other words, the WZW fixed point in (||) is unstable and (^) defines 
a massive field theory called the Gross-Neveu model. For negative g, the coupling is marginally 
irrelevant and the WZW fixed point is stable, a fact which will become important for sigma 
models with 6 = ir. 

As one would expect, there are particles in the O(N) Gross-Neveu model with the same 
quantum numbers as the fermions in (||). These are in the defining (^-dimensional vector) 
representation of O(N). For general G, this is also true: there are particles in the defining rep- 
resentation of G. However, there is more. For O(N), there are kinks as well, as follows from a 
semi-classical analysis of the action (@), These particles are in the spinor representations of O(N) 



[23]. Thus the 0(N) Gross-Neveu model exhibits charge fractionalization. For even N, there are 
two spinor representations, of dimension 2 N l 2 ~ l ; for odd N there is one spinor representation, 
of dimension 2^ N ~ 1 ^ 2 . However, there still more quasiparticles, which are bound states of the 
fermions. As will follow from studying the bound-state structure of the exact S matrix, there are 
particles in all of the fundamental representations. The fundamental representations for any Lie 
group are the representations with highest weight vectors ji 1 whose inner product with the simple 



root vectors a 3 obeys [i % ■ a? = S 11 [24]. There is one fundamental representation corresponding 
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to each node on the Dynkin diagram. The fundamental representations of SU(N) have Young 
tableau with only one column: the defining (iV-dimensional) representation, the two-index an- 
tisymmetric tensor, all the way through the N representation. Thus for G = SU(N), there are 
particles corresponding to any antisymmetric tensor representation with j indices, j = 1 . . . N — 1. 
All of these can be considered as bound states of the iV-dimensional representation. For O(N), 
the fundamental representations are the antisymmetric tensors and the spinor representation(s). 
All of these can be obtained as bound states of the spinor (s). 

The Gross- Neveu models are integrable, so their exact S matrices can be derived. For SU(N), 
the tensor product of two vector representations contains the symmetric tensor and the antisym- 
metric tensor. Thus the S matrix for two particles in the vector representation is of the form 

s vv = f vv Vs + f vv p ^ (10) 

where the projection operators on the two-index symmetric (S) and antisymmetric (A) represen- 
tations are explicitly 

v s {a i {e a )b j {Q h )) = i {oi{e a )bj[e h ) + a j {e a Me b )) 

PA(ai(Pa)bj(e b )) = \ (Oi(0 o )&#6) - o.j(9 a )bi(e b )) . (11) 

tii tii 

The subscripts in ai and bj represent the i and j particles in the vector multiplets, so the S 
matrix in (|BJ) is indeed a matrix. 

Requiring that the S matrix obey the Yang-Baxter equation means that 



f\ v _ 9 + 2mA 



(12) 



where A will be determined shortly. The simplest solution for the overall function consistent 
with crossing symmetry and unitarity for particles in the vector representation of SU (N) is 

fs V = S mYm where 

gVV _ ^ (j; ~ 2m) ^ (im + ^) (in\ 

mm r(i + ^)r(-^ + A) { j 

This is called the minimal solution, because the resulting S matrix has no poles in the physical 
strip: a zero in ( |l3| ) cancels the pole in (|l2|). It is not the unique solution, because of the CDD 
ambiguity. 

The Gross-Neveu model has bound states and so the S matrix must have poles in the physical 
strip. This means that fg = S^ l Y n X(9), where X{9) contains the poles. If there are to be bound- 
state particles with in the antisymmetric representation of SU(N), but not in the symmetric 
representation, then there must be a pole in /jf but not in f]f . The relation ( |T2] ) means 
the pole must be at 6 = 2iriA and so particles in the antisymmetric representation have mass 
mA/mv = 2 cos irA. To determine A, one must ensure that the bootstrap closes. This means that 
one does not generate an infinite number of particles: the poles in all the S matrices correspond 
to a finite number of particles. For SU(N), the bootstrap closes if A = 1/N. For general 
groups, the Gross-Neveu bootstrap closes if A = 1/h, where h is the dual Coxeter number. In 
the rest of this paper, A = 1/h, with h appropriate for the case at hand. Putting (Hj^ 11 ,12,0) 



and the bootstrap together means that the S matrix for the scattering of vector particles in the 



Gross-Neveu model is [22] 



x m *lUo)[v s + t±^v A ) .id 
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where 

S inh((fl + 27rcA)/2) 
X[e) ~ sinh((#-2™A)/2)' (15) 

The extra factor X{9) indeed reinstates the pole at 9 = 2iriA canceled by S m i n . Using this boot- 
strap to compute the S matrices of the bound states, one finds for example that the scattering 
matrix S AV for a particle in the vector representation with one in the antisymmetric represen- 
tation has a pole at 9 = 37rzA, leading to particles with mass m^jmy = sinh(3A7r)/ sinh(A7r). 
These particles are in the representation with highest weight ^3 (the Young tableau with 3 boxes 
and one column, or equivalently the three- index antisymmetric tensor). 

One keeps repeating this bootstrap procedure and finds particles in all the antisymmetric 
representations, ending up with the particles in the conjugate representation N representation 
arising from the bound state of N — 1 vector particles. This provides a substantial check on the 
S matrix, because the crossing relation (P) relates S vv to S vv . One can indeed check that the S 
matrices built up from ( |l4|) satisfy this relation. For the SU(N) Gross-Neveu model, there is one 
copy of each fundamental representation pi] . The dimension of the antisymmetric tensor with j 
indices (j boxes in the one-column Young tableau) is N\/j\(N — so there are 2^ — 2 particles 
in all. Each multiplet of particles has mass M smirj/N, where M is an overall scale related to 
the coupling constant g. The particles in the N — j representation are the antiparticles of those 
in the j representation. One can check that there are no additional particles by computing the 
energy in a background field as done in subsequent sections p3] , or by computing the free energy 
at non-zero temperature [26]. Moreover, the SU(N) Gross-Neveu model can be solved directly 
using the Bethe ansatz p2]]; the results agree with those above. 



The S matrices for the 0(N) Gross Neveu model are derived in the same manner as those 
with SU(N) symmetry, but there are a number of additional complications. 

The scattering of particles both in the vector representation of O(N) is of the form 



s v v = jrv v Vs + ^ v Va + J7VV Pqj ( 16) 

where the projection operators on the symmetric (S), antisymmetric (A) and singlet (0) repre- 
sentations are explicitly 

1 1 N 

Vs(aibj) = - (flibj + ajbi) - —Sij a k b k 



N 

k=l 



V A (aibj) = - (ciibj - ajbi) 

N 



Vo(aibj) = jjdij^Takbk (17) 
fc=i 

where I have suppressed the 6 dependence. For example, the scattering process a\b\ — > a%bi has 
S matrix element {{N -2)FY, V + NF\ V + 2F^ v )/2N, while the process a x b x -> a 2 b 2 has element 
{J~q V — J-^ V )/N . The extra term in (|i~6| ) as compared to ( |To| ) stems from the fact that the trace 
Sijdibj is an SO(N) invariant. 

Requiring that the S matrix obey the Yang-Baxter equation means that Q 

TY_ + 2mA 

T v s v 9 -2mA' { ' 
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FX V _ + 27TiA 8 + ITT 

T v s y 9 -2mA 9 -m [ ' 

The minimal solution (no poles in the physical strip) for the overall function for 0(N) with 
particles in the vector is !F^ V = S^J n , where 

cW (fi\ — E (j: ~ E (j in) E (jzH A" ^} E ~ M it ^) fofll 

mmU r(i + ^)r(i-^)r(-^ + A)r(i + ^ + A) 1 j 



The S matrix with /J y = S^J n defined by (g^ll® is the S matrix for the 0(N)/0(N - 1) 
sigma model []|]. The S 1 matrix does not contain any poles because the 0(N)/0(N — 1) model 
has no particles in any representations other than the vector, and hence no bound states. 

The O(N) Gross-Neveu model has bound states. Just like in SU(N), there are particles 
in all the antisymmetric representations, whose S matrices follow from the bootstrap. These 
representations are self-conjugate: there is no N representation. That means that both the bound- 
state (s-channel) pole at 9 = 2niA and the crossed pole (i-channel) channel at = m — 2itiA 
appear in Sq^, yielding 

- *CW* - e )&) (v s + v A + ^ Po ) (2D 

where X(9) is as in (pl|) with A = 1/h, with the dual Coxeter number here h = N — 2. This 
value of A ensures the bootstrap closes. The particles obtained by fusing the vector particle j 
times therefore have mass rrij/mv = 2 sin(j7rA)/ sin(7rA). One difference between the O(N) and 
SU(N) Gross-Neveu models is apparent in (|2~I|): if there is a pole in the antisymmetric channel 
at 9 = 2iriA, there is also one in the singlet channel. This means that there are extra bound 
states in the O(N) Gross-Neveu model. For example, there are iV(./V — l)/2 + 1 particles with 
mass ni2- N(N — l)/2 of them transform in the two-index antisymmetric representation of O(N), 
while one is a singlet under 0(N). In general, at mass m,- there are particles in all the A:-index 
completely antisymmetric representations k = 0, 2, . . . j for even j and k = 1,3, ... j for odd j 



[15]. Despite this additional complication, one can in principle obtain the S matrices for all of 
these states by using fusion. This was done explicitly for k = 2 in ]17| . 

The S matrix for the kinks in the spinor representation in the 0(N) Gross-Neveu model 
is given in j jX5| , [l6| , 27] . It can be written out explicitly, but since the answer is somewhat 
complicated I will not give it here. It is derived by using the fact that kink-kink bound states 
contain the particles in the antisymmetric representations as bound states; in fact is a sort of 
reverse bootstrap. The mass of the spinor states is given by m sp i nor /m v = 1/ sin(-7rA). 

In general, the particles of the G Gross-Neveu theory are in fundamental representations of 
G. Their masses and multiplicities are known from considering the bound-state properties of the 
exact S matrix. 



5 The particles and their scattering matrices 

In this section I give the particle content and exact scattering matrices of the SU(N)/SO(N) 
and 0(2P)/0(P) x O(P) sigma models. In the next section, I will use the exact S matrix to 
compute the energy in a background field. 

The representations of the global symmetries of sigma models are somewhat more elaborate 
than that in the Gross-Neveu model. The results are different for each model, and also change 
dramatically if certain extra terms are added to the action (S) . 



11 



5.1 SU(N)/SO(N) sigma model 

For SU(N)/SO(N), the action is written in terms of a symmetric unitary matrix field <£. Thus as 
opposed to the Gross-Neveu model, one expects there to be particles forming the symmetric rep- 
resentation of SU(N), which has dimension iV(iV + 1)/2. Indeed, it was long ago established that 
for N = 2 (the sphere sigma model), the particles are in the symmetric (triplet) representation 
of SU(2) H. For any N, there are non-trivial non-local conserved currents in this sigma model 
p8| . These conserved currents are consistent with particles in the symmetric representation, and 
also with bound states in all representations of SU(N) with Young tableau with two columns and 
rectangular (i.e. the same number of boxes in each of the columns). In group-theory language, 
these are representations with highest weight 2^ (the fundamental representations arising in the 
Gross-Neveu model have highest weight /i J ). 

Building the bound-state S matrices with the bootstrap relation ([?]) is one example of a 
general procedure called fusion [29|. Fusion exploits the fact that the S matrix can be written as 



the sum of projection operators to find new solutions of the Yang-Baxter equation. It is useful 
for things other than just the bound states in a given theory. The solution of the Yang-Baxter 
equation S ss for particles are in the symmetric representation can also be obtained from S . 
The solution is obtained from (^) like before, using instead the pole at 9 = — 2-7riA. The tensor 
product of two symmetric representations contains three representations with Young tableaux 



4 /"i = i i i n 2 a*i + /"2 



2^2 



It is convenient to label particles in the symmetric representation with two indices and the 
constraint ajj = aj%. The projection operators are then explicitly 

T^Vi (aijbki) = - (aijbki + ciikbji + aubjk + akihj + a^bn + (ijibik) 

'P2fi 2 ( a ijbki) = -(^aijbki + 2akihj — a ik bji — aiibjk — aj k bii — ajib ik ) (22) 

This means, for example, that the S matrix element for scattering an initial state of a\i and 613 
to a final state of a\2 and 613 is (2/^ + 3/£^ + + /2^)/6, while the element for scattering the 
same two particles and getting a\\ and 623 in the final state is (/£fj — /^)/6- 

This enables us to find the S matrix for the particles in the SU(N)/ SO(N) sigma model at 
9 = in the symmetric representation: 

S"» = X WS ™ (0) ( V lm + j-^ 7W + j^-k j-^ «*■ ) • < 23 > 



where X(9) is the same as in the Gross-Neveu model. The overall function <S^f n is determined 
by using unitarity, crossing and the bootstrap. The result is that 

ss 9-27riA r (1 - r + 2A) 

mmU ^ + 2 - A r(i + ^)r(-^ + 2A)- 1 j 

The symmetric representation has highest weight 2\i\. The pole in the S matrix at 9 = 2niA 
results in a bound state transforming in the representation with highest weight 2/U2, with mass 
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rn 2^ 2 / rn s = 2cos7rA. In fact, the factor X ensures that the spectrum of bound states is just 
like that in the SU(N) Gross-Neveu model. Building the S matrices for the bound states in the 
SU(N)/SO(N) sigma model proceeds just like the SU(N) Gross-Neveu model. For example, 
the prefactor in ( |23| ) ensures that S determined by the bootstrap is the same as that following 
from crossing. To close the bootstrap, A = 1/h as before, with h = N for SU(N). There are 
particles in all representations with highest weight 2/ij. Each representation appears just once 
with mass is M sinirj/N. For N = 2, this model reduces to the 0(3)/0(2) model (jl|), and the S 
matrix (^) becomes that of a massive 0(3) triplet, familiar from [g]. 

The SU(N)/SO(N) sigma model with no 6 term is therefore an integrable generalization of 
the sphere sigma model to all SU (N). I will now show that the behavior at 8 = tt also generalizes. 

To find the particles and scattering matrices when 9 = 7r, it is first necessary to understand 
the non-trivial fixed point. Generalizing from N=2, it is natural to assume that it is SU(N)i. 
In fact, the argument from [fj], ^] can be adapted to show that this is completely consistent with 
the symmetries of the problem. To apply the consistency argument here, it is first useful to study 
the symmetries of the SU(N) Gross-Neveu model. To avoid confusion with the sigma model, I 
denote the field in the Gross-Neveu model as w. With the definition (^) as a perturbed WZW 
model, w must be an SU(N) matrix. Therefore, it transforms under the SU(N) symmetry as 
w — > UwW , where U is an element of SU(N). The currents also transform as j — ► UjU^. In 
other words, w and j are in the adjoint representation of SU(N). The action (g) is manifestly 
symmetric under the symmetry SU(N)/Zin ■ The Z^r is the center of SU(N), consisting of of 
matrices £11, where is a Nth. root of unity, and I is the identity matrix. The matrices f2I 
commute with all elements of SU(N). The reason this discrete subgroup is divided out is that 
if U = £11, w and j are left invariant. Thus the symmetry acting manifestly on the action is 
SU(N)/Zn, not SU(N). However, the full symmetry of the Gross-Neveu model is larger. This 
can be seen in several ways. First of all, note that there are particles in the model in the vector 
representation. The full SU(N) group acts on the vector representation non-trivially. Indeed, one 
can assign an extra Z jy charge to the particles of the model, defined so that the particles in the 
j'-index antisymmetric representation have charge IP . Another way of seeing this extra charge is 
by examining terms which do not appear in the action (^). The £[/(iV)/Zjv-symmetric operators 
in the SU (iV)i conformal field theory are tr w 3 , where j = 1 . . . N — 1. Some of these are relevant 
operators, and if added to the action would change the physics considerably. However, they are 
forbidden from the action if the model is required to be invariant under the symmetry w — > Clw. 
This discrete symmetry is not part of the SU(N) acting on the w field, but rather makes the 
full symmetry of the model Z^r x SU(N)/Zn. This symmetry of the action and the full SU(N) 
of the particle description are completely consistent with each other if the particles are kinks in 
the w field. (They are kinks in the 1 + 1-dimensional picture, vortices in the 2 + O-dimensional 
picture.) The presence of kinks is easy to see. The field values w = I and w = £1 are not related 
by a continuous symmetry, so field configurations with w(x = —oo,t) = I and w{x = oo,t) = QI 
are topologically stable. These are the kinks. The extra Z^r symmetry is then indeed the discrete 
kink charge. This sort of symmetry should be familiar from the sine-Gordon model, where the 
soliton charge is not an explicit symmetry of the action. 

The symmetries of the SU(N)/SO(N) sigma model at = n are similar to those in the Gross- 
Neveu model. The Zjy center acts non-trivially on the SU(N)/SO(N) matrix as seen in (||). 
Thus the symmetry of the sigma model action is the full SU(N) (or to be precise, for even N it is 
Z2 x SU(N) /Z12). Therefore it is consistent for the low-energy fixed point of the SU(N)/SO(N) 
sigma model to be SU(N)i. The reason is the same as in the SU(N) Gross-Neveu model: the 
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full SU(N) symmetry of the action in terms of implies a Z^r x SU(N)/Zn symmetry of the 
low-energy effective action involving w. This forbids the relevant operators from being added 
to the SU(N)i action. This argument shows that the effective action near the low-energy fixed 
point contains only 5C/(A^)-invariant irrelevant operators. The operator jlJr is SU{N) invariant 
and of dimension 2. Thus the effective action here is like the Gross- Neveu action (|9|), except here 
the coupling g is negative so that the perturbing operator is irrelevant. The coupling g is related 
to the scale M; when M — > oo, g — > from below, and the model reaches the low-energy fixed 
point. 

The form of the effective action near the low-energy fixed point gives the S matrices Sll and 
Sun immediately, because as explained above they are independent of M and therefore follow 
solely from SU(N)\. They must be in fact those of the massless limit of the Gross-Neveu model. 
Written as functions of rapidity they are 

St L (9) = S a B b R (9) = S a G b N (9). (25) 

Thus the spectrum of these gapless particles is the same as that of the massive particles of the 
Gross-Neveu model: there are particles in any antisymmetric representation, a left-moving and 
right-moving set for each representation. Even though the particles here are gapless, the mass 
ratios of the Gross-Neveu model still appear in the definition of rapidity. For example, a right 
mover in the representation \J has energy E = mje 6 '. The charge has fractionalized, as compared 
to 9 = 0. 



For N=2, the S matrix Slr(0) as a function of rapidity is the same as Sll{9) JlO[- This S 
matrix satisfies the Yang-Baxter equation, and has the appropriate SU(2) (not SU (2) l x SU (2)r) 
symmetry. However, the situation is not quite as simple for general N, because there are poles 
in Sll and Srr in the physical strip, resulting in the bound states. These poles are forbidden 
in Slr- However, it is easy to remove these and still have a sensible S matrix. The bound- 
state S matrices S G b N are labeled with a, b = 1 . . . N — 1 corresponding to the antisymmetric 
representations with a and b indices. For example S 11 = S vv and S 21 = S AV . The unwanted 
poles in S ab all arise from the factor X ll {9) = X{9) and its fusions. The matrix Sq N contains 
the overall function 

a b 

X ab (9) = Y\ Yl X (9 + [2{i +j-l)-a- b]mA) . (26) 
i=\j=\ 

For example, X 12 = X(9 + 7riA)X(9 — iriA). The prefactor X ab contains all the poles in the 
physical strip in S^ N . It is easy to check that X ab {9)X ah {-9) = 1 and that X M (9) = X ab (m-9). 
Therefore 

St R {0) = S a M0)/X ab (9) (27) 

satisfies crossing, unitarity, the Yang-Baxter equation, and has no poles in the physical strip. 
Thus this is the S matrix for left-right scattering in the SU(N)/SO(N) sigma model at 9 = it. 
For JV = 2, X = 1, and the result reduces to that in JlO[ . 

To prove that this picture is correct, in the next section I will show that these S matrices 
give a model which is SU(N)\ in the low-energy limit, but is the SU(N)/SO(N) sigma model at 
high energy. Moreover, in the corresponding c function is calculated, and indeed flows from 
the c = (N - 1)(N + 2)/2 (the value at the trivial high-energy fixed point in SU(N)/SO(N) to 
c = N — 1 (the central charge for SU (N)i). 
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5.2 0(2P)/0(P) x O(P) sigma model 

As with the Gross-Neveu models, the 0(2P)/0(P) x O(P) sigma model is a slightly more com- 
plicated version of its SU(N) analog, the SU(N)/ SO(N) sigma model. The field $ is symmetric 
and unitary (and also real and traceless), so again one expects particles in the symmetric repre- 
sentation of 0(2P) (which is P(2P + 1) — 1 dimensional), and its various bound states. This is 
what happens for P = 2, where the model reduces to two copies of the sphere sigma model: there 
are six particles: three in each of the symmetric representations. Here, for P > 2, the non-local 
conserved currents have not been found, although some interesting results for the local currents 
in the classical model (g small) have been found |}^] . 

The S matrices for the 0(2P)/0(P) x 0{P) model can be constructed by fusing the 0(2P) 
Gross-Neveu results. The S matrix for particles in the symmetric representation of 0(2P) (highest 



weight 2/xi) is of the form [17] 



,„s 9 + 4mA 9 + in + 2ni A _ \ , . 

9 + Am A 9 + 2mA / 9 + in + 2mA 9 + in + 2mA 9 + in 

+ 9-4niA 9 -2niA \ 2M2 + 9 - in - 2mA M2 + 9 - in - 2niA 9 ^ ' " 



in 

where A = 1/h as always, with h = 2P — 2 here. The minimal solution (no poles in the physical 
strip) for the prefactor is 

ss 9 - 2mA r(i-^)rQ + ^)r(^ + 2A)r(i-^ + 2A) 

mmU e + 2 - A r(i + ^)r(i-^)r(-^ + 2A)r(i + ^- + 2A)- { j 

Like the 0(2P) Gross-Neveu model, to get bound states there must be poles at 9 = 2niA and 
9 = ni{\ — 2 A). This means that the prefactor of S ss is 

R{9) = X(9)X(in-9)S^ m {9). (30) 

Because the factor X(9)X(in — 9) is the same, the bootstrap for the 0(2P)/0{P) x O(P) model 
gives bound states with the same spectrum as the 0{2P) Gross-Neveu model. Also like the 
0(2P) Gross-Neveu model, there are particles which do not follow obviously from the action. 
These are in "double-spinor" representations, formed by the symmetric product of two spinor 
representations. In group theory language they have highest weight 2/i s , where fi s is the highest- 
weight of a spinor representation; they are {2P — 1)!/P!(P — 1)! dimensional. These particles 
presumably are kinks like in the Gross-Neveu model, but it is not clear how to extract this 
information directly from the action. The S matrix for these kinks is quite complicated, since 
many representations appear in the tensor product of two double-spinor representations. It can 
presumably be obtained by fusion of the spinor S matrices, or by the reverse bootstrap like the 
Gross-Neveu kinks. The presence of these particles is confirmed by studying the free energy at 



non-zero temperature [31]. 



Therefore, in both series of sigma models there are particles in all representations with highest 
weight 2//. However, here there are even more degeneracies and more representations appearing. 
Since there is a pole in ( [29] ) at 9 = 2niA, then there must be bound states not only in the 
representation with highest weight 2/i2, but also in the singlet and in the antisymmetric rep- 
resentation (highest weight ^2)- Thus charge is fractionalized even at 9 = 0. For example, in 
0(8), the vector and spinor representations are 8-dimensional, the antisymmetric representation 
has dimension 28, the symmetric and double-spinors dimension 35, and the representation with 
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highest weight 2^2 has dimension 300. Thus in the 0(8) Gross-Neveu model, there are 8 particles 
in the vector (mass 2Msin-7r/6 = M) 8 particles in each of the spinor representations (mass M), 
and 28 + 1 particles of mass 2Msin7r/3 = y/3M. Thus there are 53 stable particles in the 0(8) 
Gross-Neveu model. In the 0(8)/0(4) x 0(4) sigma model, there are 35 particles in each of the 
vector and two double-spinor representations, all of mass M. There are 300 + 28 + 1 particles of 
mass y/3M, giving 434 particles in all. 



The behavior when 9 = n in the 0{2P)/0{P) x O(P) model is also reminiscent of the 
SU(N)/SO(N) model at 9 = it. The arguments follow in the same fashion. The model has the 
0(2P)\ WZW model as a stable low-energy fixed point. The 0(2P)\ WZW model is equivalent 
to 2P free Majorana fermions, or equivalently 2P decoupled Ising models. The word "free" 
is slightly deceptive, because just as in a single 2d Ising model, one can study correlators of 
the magnetization or "twist" operator, which are highly non-trivial. The consistency argument 
is simpler here: the only relevant 0{2P) symmetric operators are the fermion mass and the 
magnetization operator; neither is invariant under the symmetry <3? — > — of the sigma model. 



The c- function is computed in [31], and indeed flows from c = P 2 to c = P as it must. 

The S matrix for the quasiparticles follows from the Gross-Neveu model. The same arguments 
applied above to the SU(N)/ SO(N) sigma model at 9 = tt show here that the massless left- 
and right-moving particles have the same spectrum as the 0{2P) Gross-Neveu model. Charges 
fractionalize: the left- and right-moving particles are in the fundamental representations of 0{2P). 
The S matrices for the particles in antisymmetric representations are 

st L {e)=s<$ R {e) = s£ N {6) (3i) 
s a Me) = s$ N {o)l{x a \e)x ab {m - e)) (32) 



For particles in the spinor representations, the poles can easily be removed as well; sec [16] for the 
definition of X ab for the spinor representations. Even though the field theory at the low-energy 
fixed point is free fermions, the S matrix is factorizable away from the fixed point only in the 
basis related to the Gross-Neveu model. 



6 Matching perturbative expansions 

In this section I compute the energy of these sigma models at zero temperature in a background 
magnetic field. This is very useful for several reasons. First of all, it allows a direct comparison of 
the S matrix to perturbative results. This in particular ensures that the S matrices are correct as 
written in the last section. For example, it eliminates the possibility of extra CDD factors and/or 
extra bound states. Second, because the effect of a 9 term is non-perturbative, the perturbative 
expansions for the sigma models at 9 = and tt must be the same. Thus even though the S 
matrices for 9 = and tt are very different, the energy must have the same perturbative expansion. 
I verify that this is true for the above S matrices. 

The abelian subgroup of the group G is U(l) r , where r is the rank of the group. Thus a 
model with a global symmetry G has r conserved charges. These charges can be coupled to 
a background field, which is constant in spacetime. In the sigma models 0(2P)/0(P) x O(P) 
and SU(N)/0(N) where $ is a symmetric matrix, this means that the Euclidean action @ is 
modified to 

S=-tr I d 2 x (V$ f - A T & - $ f i T ) + i$ + $i) (33) 
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where A is a matrix in the Cartan subalgebra of G (the Cartan subalgebra is comprised of 
the generators of the abelian subgroup of G). For SU(N)/SO(N), A is diagonal with entries 
(A 1 ,A 2 ,...A N ) and the constraint Y.i A % = °- For 0(2P)/0(P) x O(P), the matrix can be 
written in the form A jk = Yd=iAi^j,2l-lh,2l ~ S j,2l S k,2l-i)- 

Because A has dimensions of mass, the energy depends on the dimensionless parameters 
Ai/M or Ai/M. The strength of the background field controls the position of the theory on its 
renormalization group trajectory and, in particular, in the limit of large field the theory is driven 
to the ultraviolet fixed point. To be near the UV fixed point, all non-zero Ai/M must be large 
in magnitude, so let A be one of the non-zero components. The perturbative expansion around 
this fixed point therefore is an expansion for large A/M. The computation of this expansion is 
a fairly standard exercise in Feynman diagrams; for details closely related to the cases at hand, 
see [25]. The only effect of the sigma model interactions resulting from the non-linear constraints 
(|3|) to the one-loop energy comes through the running of the coupling constant. Specifically, the 
two-loop beta function is of the form 

P(g) = = -pVM - pVM - • • • 

where (3\ and 02 are model-dependent, but universal for these sigma models. Solving this equation 
means that 

-j-r = ft In A/ A + ^ ln(ln A/ A)) + 0(1/ ln(A)) (34) 
9{ A ) Pi 

where the scale A depends on the perturbative scheme used. The zero-temperature energy for 
SU(N)/SO(N) through one loop is 

E{A) - £7(0) =--rfi E(^') 2 " h ~ ( ln |A ~ Ajl 

with a similar formula for 0(2P)/0(P) x 0(P). For both cases, this means that at large A/M, 
the energy is of the form 



E(A) - E(0) oc A 2 {ln(A/M) + \n(\n{A/M)) + . . .) (35) 

(Pi) 



In particular, note that the ratio of the first two terms is universal. These logarithms are char- 
acteristic of an asymptotically free theory, where the perturbation of the UV fixed point is 
marginally relevant. 

I now explain how to calculate the ground-state energy E(A) — E(0) directly from the S- 
matrix, giving the promised check. In this picture, the model is treated as a 1 + 1 dimensional 
particle theory at zero temperature. Turning on the background field has the effect of changing 
the one-dimensional quantum ground state. If one were working in field theory or in a lattice 
model, this would change the Dirac or Fermi sea. In the exact S-matrix description, a similar 
thing happens: the ground state no longer is the empty state — it has a sea of real particles. 
For example, a particle aj in the vector representation of SU (N) has its energy shifted by Ai in 
a magnetic field. If the total energy is negative, then it is possible for such a particle to appear 
in the ground state. Because the technicalities are slightly different, I treat the cases of massive 
and massless particles separately. 
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6.1 Massive particles 

For simplicity, I first treat only the case where the magnetic field is chosen so that only one kind 
of massive particle of mass m and charge q appears in the ground state. In integrable models 
obeying the Yang-Baxter equation, the particles fill levels like fermions: only one can occupy a 
given level. Then the ground state is made up of particles with rapidities 

-B <0 1 <9 2 < ...9m < B 

for some maximum rapidity B. When the length L of the system is large, these levels are very 
close together, so the density of particles p{6) is defined so that p(9)d9 is the number of particles 
with rapidities between 9 and 9 + d9. The energy density of the ground state is then 

1 r B 

E{A) -E(0) = -V / d9 p(9) (- A + m cosh 9). (36) 

L a J - B 

I have normalized A so that the charge q = 1. The exact 5-matrix allows us to derive equations 
for the p{9) and B. Imposing periodic boundary conditions on the box of length L requires that 
the rapidities 9i of the particles in the ground state all satisfy the quantization conditions 



Y[S{0 i -6 j ) = l (37) 



for all i. This is an inter acting-model generalization of the one-particle relation m-sinn^ = Pk = 
2irrik/L to the case where the particles in the ground state scatter elastically from each other 
with S matrix element S(9). In the large L limit, taking the derivative of the log of ( |37| ) gives 
an integral equation for the density p(9): 



mL 



2vr 

where 



B 

p(9) = ^ cosh 9 + / d9' p(9')(j)(9 - 9'), (38) 

B 



^ K ' 2vr 89 

This equation is valid for \9\ < B; for \9\ > B, p{9) = 0. The maximum rapidity B is determined 
by minimizing the energy equation (|3£|) with respect to B subject to the constraint (^). If the 
particles were non-interacting, the density would be mL cosh 9 and m cosh B = A, but the effect 
of the interactions is quite substantial. 

These equations can be put in a convenient form by defining the "dressed" particle energies 
e(0) as 

f B 

e{9) = A-m cosh 9+1 d9'(f)(9 - 9')e(9'). (39) 

J-B 



Substituting this into (pq) and using (38 yields 

E(A) - E(0) = / d6 coshfl e{9). (40) 

In this formulation, B is a function of A/m determined by the boundary condition e(±B) = 0. 

Consider the 0(2P)/0(P) x O(P) sigma model, and choose the magnetic field to be A\ = A, 
Ai = for i 1. The massive particles used above (a,ij = ciji) are not eigenstates of the 
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matrix magnetic field A, but one can easily change basis to particles which are. This merely 
requires taking linear combinations of particles of the same mass. The state which has the 
largest eigenvalue of this magnetic field is d = an — 022 + 2ia\2- This basis also has the advantage 
that scattering of d particles amongst themselves is diagonal (as opposed to the scattering of 
an with bn, which does not necessarily yield an and bn in the final state). The S matrix 
element for scattering d(9i)d(02) to d{6\)d{92) is R(Qi — 9 2 ) as given in (p0|). The state d has the 
largest eigenvalue of the magnetic field, and is the only state with this eigenvalue. I therefore 
make the assumption that it is the only kind of particle which appears in the ground state. This 
assumption is standard in these computations; it can be justified by a careful treatment of the 
zero-temperature limit of the finite-temperature density equations discussed in |3lJ]. The kernel 



of the integral equation (39) is then proportional to the derivative of log R(9). It is convenient 



to write this in Fourier space as 

<j)(0) = I duo e lwU {l - K{uj)) (41) 



OC 



_„ AM 2 cosh((l - 2A W2) sinh(27rAH) 

M w ) = e w 77^ ( 42 i 

cosh(7ru;/z) 

where A is as always 1/h, with the dual Coxeter number h = 2P — 2 for 0(2P). 

For the SU(N)/SO(N) sigma model, the states ay are eigenstates of the magnetic field 
operator, with eigenvalue Aj_ + Aj . If we chose the magnetic field to be A\ = A, Aj = —A/ (N — 1) 
for j > 1, then the particle an has maximum charge. Again making the assumption that this is 
the only particle in the ground state, the Fourier transform of the resulting kernel is 

2c -7rA\u,\ sin H( l ~ A)vr|tj|)sinh(27rAa;) 



sinh(7ra;) 

However, it is convenient here to instead choose the fields to be Ai = —A2 = A, Aj = for j > 2 
(this choice of field was useful also in the supersymmetric CP n model [p3| ) . Then there are two 
particles with largest eigenvalue: an and the antiparticle (in the iV representation) 022- Thus 
they both appear in the ground state. These two particles scatter diagonally among themselves 
and each other. The S matrix element Si for onon — > anau and 022^22 — ► ^22^22 is 

s 1 (e) = s%& l (o)x(p). 

The element S2 for 011022 - ► 011022 is found easily from (^) by using crossing. It is 



(in - - 2n)(iir - - /i) 



Because of the symmetry between the two kinds of particles, their ground-state densities must 
be the same. The ground-state energy then follows from a simple generalization of the above 
analysis: the equations (39) and (EH) still apply, with 



here. Plugging in the explicit expressions of the functions yields the useful fact that S\(9)S2(9) = 
R(9), where R(9) is the function appearing in ([?(]). Thus with these choices of magnetic field, the 
0(2P)/0(P) x 0(P) and the SU(N)/SO(N) sigma models can be treated by using the kernel 
the only difference is the A = 1/(2P — 2) in the former and A = 1/N in the latter. 
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The linear integral equation (39) cannot be solved in closed form. However, there is a gen- 
eralized Weiner-Hopf technique which allows the perturbative (and non-perturbative) expansion 
for large A/M to be obtained ]3q| . I discuss this technique in the next subsection. In order to 
compare the energy in the 9 = sigma model with the perturbative computation (|3~5|), I can rely 



on the results of 25]. There the first few terms in the large A expansion for models with kernels 
like ( f42| ) are given. The technique requires that the kernel fC(u>) be factorized into 



1 



(43) 



where K. + {uj) has no poles or zeroes (and is bounded) in the upper half plane Im{oj) > 0, while 
JC-(uj) = /C+(— u>) has no poles or zeroes and is bounded in the upper half plane. Then 



where 



2iAuj 



TT 



r(|-i|) 



(44) 



2A ln(2A) + (i - A) ln(i - A) + i In 2 - A. 



The factor e 1 ^^ ensures that JC+ is bounded appropriately as \u>\ 
When for small £, K + {i^) goes as 



oo in the upper half plane. 



k 



1 + ...) 



the energy at large A is 25] 



E(H)-E(0) = --A 2 



+ 



This expansion is of the same form as (|3^). The two are equal if 



02 



{Pi? 



The explicit kernel 



yields 



s + 



A. 



1 



(45) 



(46) 



while perturbative computations (see e.g. [34] and references within) give for SU(N)/SO(N) 

2 N + N 2 /2 



{0if 

while for 0(2P)/0(P) x 0{P) the ratio is 

02 



N 2 



2P 2 - 2P 



(/3i) 2 (2P-2) 2 

Using the relations A = 1/N and A = 1/(2P — 2) respectively, one indeed sees that that the 
condition (|4^) holds. This is a substantial check, and coupled with the fact that S matrix also 
gives the correct central charge, leads me to be completely convinced that the 5 matrices discussed 
above are indeed the correct sigma model S matrices. 
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6.2 Comparing the perturbative expansions at 6 = and n 

In this subsection I analyze the equations for the ground-state energy in a magnetic field in more 
detail. The final result will be that the entire large- A perturbative expansions for 8 = and for 
8 = ir are identical, but that non-perturbative contributions differ. This effectively confirms the 
S matrices and spectrum above, and the identity of the low-energy fixed points. 

The generalized Weiner-Hopf technique is discussed in detail in the appendix of |pq| . The 
equation to be solved is of the form 

e(9)- [%(8-8')e(8')=g(8). 

J-B 

where e{9) and g{9) are both vanishing for \9\ > B, and (f)(8) = <p(—8). Defining the Fourier 
transforms g(w) and K(uj) (the latter related to 4> as in (|4l|)) gives 

oo 

ILO0 



e-^{?H/C( W )-M = 0. 

This equation is valid for \9\ < B. Fourier-transforming this gives 

e(w)/C(w) - g(u) = X + (uj)e lujB + X_e^ B . (47) 

The functions X± arise, roughly speaking, from the analytic continuation of e(8) to \9\ > B. The 
extra factors e lujB ensure that X+(uj) is analytic in the upper half plane, while X—(u) = X + (— uj) 
is analytic in the lower half plane. The relation ( |47| ) can be split into two equations, one involving 
poles in the upper half plane and the other the lower. To split a given function, one uses 

/(w) = + [/(«)]_ 

where 



2-7rz J-oo u' — uj =F i5 

where 5 is a positive real number tending to zero. The functions [f]± are analytic in the upper 
(+) and lower (— ) half planes. Because e(9) and g(8) are zero for \8\ > B, the functions 

e±(o;) =e(u)e^ B 9± {l>) = g(uj)e^ B 

are similarly analytic in the upper half and lower half planes. Equations for the functions X± 
can be derived by exploiting these analyticity properties and the factorization (|43|). Namely, (|47| ) 
implies the equations 

= [ 9± M/C T M]± + [X ± (cu)/C T (c)e ±2 - B ] ± , (48) 

while the X± are given by 

X ± {u)K±{u) + [X T (cj)IC ± (u;)e^ B ] ± = -[g T (u)IC ± (u;)] ± (49) 

Unfortunately, it is not possible to solve the equations for X± explicitly. However, this form does 
allow the l&rge-A expansion to be systematically developed. 
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To study the large A expansion, it turns out to be much easier to consider the more general 
kernel 

cosn (T 7rAw ) 2cosh((l - 2A)7rt i ;/2)sinh(2A7Rj) 

K.[ui) — — — — — - — — — . (50) 

sinh((7 + l)7rAa;) cosh(7ru;/2) 

In the limit 7 — > 00, K goes back to the sigma model kernel K. In the N=2 case, this deformation 
corresponds to deforming the sphere sigma model into the sausage sigma model |Tl]]. This kernel 
factorizes as K = 1/(K + K^), giving 

+( ^ = V 7T1 r H(7 + i)A.)r(i-, f ) • (51) 

where 

v = \i + 7A 111(7/(7 + 1)) - A In A + A. 



The equations ( |49| ) can be written as a single one by exploiting the relation X+(u) = X_(— lj). 
Defining 

K + (oj) uj + id 2 lo — i 

they become 

. . iAK+(0) ime B K+(i) f e 2lJjiB . , ,^duj' 

v{uj) = 2J V Z + ^ 111 + ; -a(J)v(J)— 52 

v ' lo + i5 2 u-i J c+ u + u' + iS y ' v J 2m y ' 

where a(u) is defined as 

, , K-(lo) 



K+(oj) 

and the integration contour circles all singularities on the positive imaginary axis. The boundary 
condition e(±B) = becomes 

iM_( )-— K.(-i) = I e 2ujW a(u;)v(Lo)— (53) 



c+ 



2iri 



while the energy (E0) becomes 



E(A) - E(0) = K 



AK + (Q)-—K + (i)l c —a^M^ 



(54) 



The contour here includes the double pole at u = i, and the poles in a(u). The function v(u) 
has no poles in the upper half plane except the explicit one at u = i. 

These equations are convenient for deriving the expansion for large A. The reason is that 
when A is large, the range of rapidities allowed in the ground state is large, so B is large (recall 



that if the particles were free, A = m cosh B). In this limit the integrals in (52,53,|jJ) are small 
corrections, and their effect can be treated iteratively. This iterative expansion is worked out 
in detail in |pql . The first correction to v(u>) comes from approximating the integral using the 
leading pieces of v{oj) in the integrand. The boundary condition ( |53| ) relates B to A. Thus the 
contributions to this integral come from the pole in v(u) at u = i and the poles in a{uj). The 
poles of a are at u = i(2n — 1)/(2 7 A), u = {2n — T)h/(2A + 1) and u = m/(2A), for n a positive 
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integer. Using the iterative procedure it is straightforward to see the form of the expansion for 
large A: it is a power series of the form 




The first series is a result of the first series of poles in a; the coefficients aij depend on the residues 
at these poles. The order A 2 corrections are a result of the pole at u = i and the other poles 
in a. Recall that 7 is large, so these contributions are much smaller. The double pole at u = % 
results in an ^4-independent piece, which is identified as E(0) 

The sigma model result of interest is recovered in the limit 7 — > 00. This complicates matters 
considerably for the power series in (|5^), since the exponent is going to zero. This is precisely 



what happens for example in the anisotropic Kondo problem as the anisotropy is tuned away [41] 
What happens in this limit is that the power series turns into an expansion in ln(^4), ln(ln(^4)), 



Note from (34) that 1/g oc ln(^4) + . . ., so the series depends on powers of g, while the order A 2 
corrections depend on e - const /9_ Thus the first series consists of the perturbative contributions to 
the free energy, while the second part is non-perturbative: the latter will never be seen in standard 
sigma model perturbation theory in g. The exact form of the perturbative expansion cannot be 
displayed in closed form, but one can build it piece by piece; the first pieces are displayed in (EET) . 



The analogous equations for E{A) for the sigma models with 9 = tt follow from their S 
matrices. I will show that the large A/M expansion is of the form (|55|), with the identical a n . 
This effectively proves that these S matrices are those of the sigma model at 6 = tt. The flow to 
the WZW model then follows immediately, because the S matrix manifestly becomes that of the 
appropriate WZW model in the large-M limit. 

The integral equations for massless particles in a background field are similar to the ones for 
massive particles. The massless particles are left and right moving, and because they are gapless, 
they begin filling the sea for arbitrarily field A. Therefore, right moving particles with rapidities 
—00 < 9 < B and left moving particles with rapidities —B < 9 < 00 fill the sea. For one species 
of right mover, and one species of left mover, the analysis at the beginning of this section can be 
repeated to derive the equations 

M f' B f'°° 
e R (9) = A-—e e + <pi(9 - 9')eR(9')d9' + / <p 2 {9 - 9')e L {9')d9' (56) 

1 J-00 J-B 

e L {9) = A-—e~ e + 4> 2 {9 - 9')e R {9')d9' +/ fatf - 9')e L (9')de' (57) 



2 



B 



valid for —oo<9<B and —B<9<oo respectively, with the boundary conditions 

e R (B) = e L (-B) = 0. 

The kernels are defined as 
The energy is 

E^(A) - fiW(O) = -— / e 9 e R (9)d9' (58) 

27T J-00 
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where I have used the symmetry 6l(9) = cr(— #)■ 

The two equations (56,57]) can be made into one by exploiting the left-right symmetry. In 



terms of Fourier transforms, €r(uj) = uj), giving 

e R (u)h(u) -e R {-u)H") -9R{u) = Y + (uj)e iu)B 

where the symmetry &2(— oS) = ^(w) is used, and gn{9) = A — me® /2. Note only one term is 
need on the right-hand side, because the integral in 9 space runs all the way to — oo. Getting rid 
of the €r(— u) gives 

«m ( *,(») - ) - ««m - - y + (-y- B + ^-c-) e -- B |^ (59) 



where as always = cj). This is now an equation of the form (f47|), and the same 

generalized Weiner-Hopf analysis can be applied. Defining 

k{UJ) — ki{LU) 



ki(u) 

and factorizing it as k(u) = 1/ (k+(u)k—(td)) in the usual way, one finds an equation just like 
(|52|), namely 

iAk+M hjO) ime B K + (i) hji) f f» iB , du/ 

v{uj) = -^Tfcl + + J c+ ^7T^ P{uJ H " m 

where 

k+{U) K2\0J) 

The boundary conditions and energy follows with the same substitutions. 

Now I can show that the sigma models have the same perturbative expansions at 9 = 
and 7r. Using the magnetic field described in the last subsection, the kernels 4>i and 02 follow 
simply from the Gross-Neveu S matrix and the relations (25 ,E7 . 31 .32). As before, it is assumed 



that only particles with largest eigenvalue of this magnetic field occupy this state. For the 
0{2P)/0{P) x O(P) models at = tt, there is only one kind of particle (left and right moving) 
in the ground state. The Fourier transforms of the kernels are 

~ sinh(( 7 + lfrAo;) cosh((l - 2A)ttu/2) 

Ki(iO) = 1 — 01 (UJ = — t \ T\ Dl 

v ' Y y ' sinh(77rA) cosh(7ru;/2) v ' 

~ sinh(( 7 - l)7rAu;) cosh((l - 2A)vr^/2) 

ko(U) = fflo w = —. n t~< 1—-. 62 

y ' r v ' sinh(77rA) cosh(7ra;/2) v ; 

where I have again included an extra parameter 7 to simplify the analysis. The sigma model 
kernels are recovered in the limit 7 — * 00, giving the exponential factors e 7rA ' w ' and e~ 7rA ' w ' 
respectively. Just like the massive case, the SU(N)/0(N) model has two particles in the ground 
state, and ends up with the same kernels. 

From these explicit forms, one finds remarkably enough that k(uS) in the massless case is 
identical to K(oS) in (|50| ) in the massive case. The only difference between the equations and 
and ( |60| ) is the extra function 

k\(uj) sinh((7 + l)A7ro;) 



&2(k>) sinh((7 — l)A7nj) 



(63) 
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in all three terms. This indeed means that the energy at 9 = tt is not the same as 9 = 0. However, 
this extra piece has no effect on the perturbative contributions, so the energy is still given by 
(|55|), This is because the extra piece introduces no new poles in the integrand in (|60|) (the poles in 
( p3| ) are canceled by zeros in k-/k+). Moreover, the residues at the poles at u = i(2n + 1)/(2A^/) 
are the same in the massive and massless cases (up to an overall sign), because 

fci(z(2n + l)/(2A 7 )) 
fc 2 (i(2n + l)/(2A 7 )) 

These residues are what determine the coefficients a n in (^), so indeed the perturbative expan- 
sions at 9 = and 9 = ir are completely identical. This was established for the sphere sigma 
model in p] ; here I have extended this proof to two infinite hierarchies of models. 

This completes the identification of the massless S matrices with the sigma models at 9 = tt. 
They give the identical perturbative expansion as the sigma model at 9 = 0, but the non- 
perturbative pieces differ. As an additional check, I have also computed the c-function by com- 
puting the free energy with no magnetic field but at non-zero temperature [31]. This gives the 
correct behavior, thus completely confirming this identification. 



7 Flows when 6 = tt 

The main point of this paper is that the stable fixed point in the sphere sigma model at 9 = tt 
is not an isolated instance. There are at least two infinite hierarchies of models which have this 
behavior. Thus current approaches to the problem of disordered electrons in two dimensions 



in both the replica approach and the supersymmetric approach (see || and [37] and references 
therein) are on sound footing. 

One obvious question is such flows happen in other sigma models. In this section I will show 
how to obtain the Sp(2N)/U(N) sigma model by perturbing the 0(4N)/0(2N) x 0(2N) model. 
I will use this to make a conjecture that at 9 = tt the former has the Sp(2N)\ WZW model as 
its low-energy fixed point. 

First, I will establish a flow between the two integrable hierarchies. The field in the 
0(2N)/0(N) x O(N) sigma model is an 2N x 2N traceless real symmetric unitary matrix. It is 
simple to see that the field configurations of the SU(N)/SO(N) models are a subspace of these. 
The configurations of the former can be written as 



A B 
B T D 



T nJ (64) 



where A, B, and D are real N x N matrices. A and D are symmetric, and trj4+ tvD = 0. To 
ensure unitarity the matrices must satisfy 

AA + B T B = DD + B T B = I 
AB T + BD = 

where / is the N x N identity. The subspace U(N) /SO(N) is obtained by requiring that A = —D 
and B = B T . The matrix A + iB is indeed a symmetric unitary N x N matrix, as can be verified 
from the preceding matrix relations. To get the SU(N)/SO(N) sigma model, one must in 
addition require det(A + iB) = 1. Thus one can flow from the 0(2N)/0(N) x O(N) sigma model 
to the SU(N)/SO(N) model adding the potential 

A J (tr [(B - B T f + {A + Df] + \det{A + iB)\ 2 ) (65) 
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Figure 1: The phase diagram of the 0(2N)/0(N) x 0(N) sigma model at 9 = tt perturbed by 
coupling A as in (|65|). The SU(N)i fixed point on the right is at A = oo. 

and making A large. If 6 = tt in the former model, then 8 = it in the latter. This flow is illustrated 
in Fig. 1. 

This flow between the low-energy fixed points can also be seen explicitly. This is the flow 
along the x-axis in figure 1. The 0{2N)\ fixed point is equivalent to 2N free massless Majorana 
(real) fermions. The right-moving fermions are denoted ip^, while the left movers are denoted 

where a = 1 . . . 2N . The 0{2N) symmetry currents are then 

and likewise for j^f . Because of Fermi statistics, j a @ = -f a , so there are N(N — l)/2 different 
right-moving currents. (I suppress the L,R subscripts: equations without them are meant to 
apply to both L and R.) These currents generate the 0{2N) symmetry, which at the critical 
point is enhanced to a chiral 0(2N)l x 0(2N)r symmetry. The theory of 2N free Majorana 
fermions is of course equivalent to a theory of N Dirac fermions defined by 

for a = 1 . . . N. From the Dirac fermions, one can form the U(N) symmetry currents 

■ab = -$ a ^b 

These N 2 right-moving and iV 2 left-moving currents generate a U(N)l x U (N)ji subgroup of the 
full symmetry. 

In the terms of WZW models, the equivalence between N Dirac fermions and 2N Majorana 
fermions is written as 

0(2iV)i = SU(N) l U{1). (66) 
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The U(l) symmetry in U(N) is generated by 



N 
a=l 

This can be split off because the generator j° commutes with all the generators j ab , so the 
corresponding WZW models are independent. For N = 2, this splitting is the simplest example 
of spin-charge separation: the C7(l) is the charge mode, while the SU(2)\ are the spin modes. 
The equivalence (|66|) is a simple example of what is called a conformal embedding [38]. Yet 



another name often used is non-abelian bosonization: the free fermions of 0(2N)i are written in 
terms of the bosonic fields of the WZW model. The simplest example of a conformal embedding 
is SU(2)i = U(l). This means that the SU{2)\ can be written in terms of a single boson <f>: the 
SU(2) currents j x ,j y and j z are cos (f>, sin (j> and d(p respectively. 

With the conformal embedding/spin-charge separation/non-abelian bosonization (|66[), it is 
now easy to see how to flow from 0{2N)\ to SU(N)\. One needs to add a perturbation which 
gives a gap to the U(l) charge mode but leaves the SU(N)i untouched. If we define the U{1) 
charge boson <p by j = dip, the perturbation can be written as 



S = S (2N) 1 + A J cosOl + <f> R ), 



Since the SU(N)i currents commute with j°, they are untouched by this perturbation. Therefore, 
in the limit A — > oo, the perturbed theory reaches the SU(N)i model. Thus the flow between 
the two sigma models can be seen in the WZW models as well: essentially all one needs to do is 
break the original symmetry appropriately. 

This is only one example of a conformal embedding. A complete list of all such embeddings 
is given in p9|. An embedding of particular interest here is 



SO(4A0i = Sp(2N) 1 SU(2) N (67) 

To realize this embedding explicitly, it is convenient to write the SO (4A r ) currents in terms of 
real antisymmetric AN x 4N matrices T?P where 

Explicitly, Tff = d£$f - 6?6%. Then the SU(2) subalgebra is given by the matrices 

T z = io y ®I 
T x = a x ®Z 
T y = a z ®Z 

where I is the 2A^ x 2N identity, the a a are the Pauli matrices, and Z is a fixed real antisymmetric 
2N x 2N matrix which obeys Z 2 = I. Since T x , T y and T z are AN x AN real antisymmetric 
matrices, they are linear combinations of the T a ^, so they do indeed form an SU(2) subalgebra 
of 0(AN). The subalgebra which commutes with the SU(2) subalgebra consists of matrices of 
the form 

ia y ®D 
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where X is the 2x2 identity and C and D are respectively antisymmetric and symmetric real 
2N x 2N matrices. For these to commute with T x , T y and T z , C and D must satisfy 

[C,Z}=0 {D,Z} = 

The subalgebra consisting of all matrices C and D obeying these constraints is precisely Sp(2N): 
by exponentiating C and D into some matrix P = e c+D , one finds P T Z P = Z, which is the 
defining relation of the group Sp(2N). Note that because C is real, P is not necessarily unitary. 
A unitary matrix can be obtained by exponentiating iC and D, giving the embedding of Sp(2N) 
in SU(2N). 

Therefore, in the embedding (|67[), the currents J a = iip (T^t-iip 1 (a = x,y or z) form the 
SU(2)n current algebra, while the currents %l) k {I®C)kiip l and ip k (ia y <g> D)^ 1 form the Sp(2N)i 
current algebra. Thus a flow from the SO(4N)i fixed point to the Sp(2N)\ fixed point occurs if 
one gives a gap to the SU(2)n modes. This can be done with the perturbation 



s = s 0(47V)l + a J tr {r L r R + jpi + j z l j z r ) 



When A goes to oo, the flow reaches the Sp{2N)\ fixed point. As a tangential comment, note that 
one can obtain a fermionic realization of the Sp(2N) Gross- Neveu model by perturbing 0(4N)% 
by the Sp(2N)i currents instead of the SU(2)n ones. 

The question now is if this flow implies a low-energy fixed point in any sigma model at 
= 7T. My conjecture is that it does. Namely, consider the Sp(2N)/U(N) sigma model. Field 
configurations in this coset space can be realized as a subspace of the 0(4N)/0(2N) x 0(2N) 
configurations. Matrices in the latter are of the form (|64l), where A, B and D are now symmetric 
2N x 2N matrices. To get Sp(2N)/U(N) requires the restrictions A = —D and B = B T as 
before, plus the additional restriction 

VZV = Z, (68) 

where V = A + iB. In other words, the Sp{2N)/U{N) subspace of 0(4N)/0(2N) x 0(2A^) 
consists of matrices of the form 

If V + V* i(V -V*) 
2\i(V-V*) -V-V* 

where V is a 2N x 2N symmetric unitary matrix obeying the condition (|68|). One can therefore 
obtain this model by a perturbation like (|65|). Here, the perturbation breaks the 0(4N) global 
symmetry down to Sp(2N). It is giving a large gap to the modes outside the Sp(2N) subgroup, 
effectively removing them from the theory. This is just how the flow described above goes from 
the SO (4A r )i WZW model to the Sp(2N)i model. It is thus very plausible that the general flows 
look like those in Fig. 2. 

This is why I have conjectured that the low-energy fixed point of the Sp(2N)/U(N) sigma 
model at 6 = it is the Sp(2N)i WZW model. This is certainly true at N = 1, because Sp(2) = 
SU(2) and this becomes the flow in the sphere sigma model. For N > 1, any non-trivial fixed 
point in Sp(2N)/U(N) should be a perturbation of the 9 = vr fixed point of 0(4N)/0(2N) x 
0(2N), namely 0(4N)i, which has central charge 2N. If there is a non-trivial critical point of 
Sp(2N)/U(N) when 9 = tt, it must have central charge less than 2N, which leaves only Sp(2N)i 

Further evidence for this conjecture comes from the replica limit N — > 0. The resulting model 
describes the spin quantum Hall effect |4Qj, and can also be studied in a supersymmetric formu- 



lation. It is found in [42] that certain correlators are equivalent to those in classical percolation. 
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Figure 2: The conjectured phase diagram of the 0{AN)/0{2N) x 0(2N) sigma model at 9 = tt 
with symmetry broken to Sp(2N). 



From this one can extract the density of states exponent in the disordered model, and indeed it 
agrees with the result from taking N — > for the corresponding exponent in Sp(2N)\ [43] (see 
also [44 1). There are many potential subtleties in taking the replica limit, so this is hardly a 
proof, but I take it as a good piece of evidence that the conjecture is true. 



8 Conclusions 

I have shown that several hierarchies of sigma models flow to a non-trivial low-energy fixed point 
when 8 = tt. This provides SU(N) and 0(2iV)-symmetric generalizations of the SU(2) results of 
HI, [lOl]. In another paper, I will calculate the c-functions for these models 31]. 



This result is a useful check on the currently-popular approach to disordered models in two 
dimensions, Most recent study has been based on the assumption that that the picture of [|| is 
fairly general, and this paper shows that the picture holds for the SU(N)/0(N) models (class 
CTI) and the 0(2N)/0(N) x 0(N) models (the GSE class). All the models discussed in this 
paper have WZW models as their non-trivial fixed points. However, if the model believed to apply 
to the integer quantum Hall plateau phase transition (the N — > limit of U(2N)/U(N) x U(N) 
model at 9 = tt M) has a fixed point, it does not seem likely that it is of WZW type, as argued 



in 3C[. Thus there is still a great deal of interesting physics yet to be uncovered in sigma models 
with 9 = tt. 



My work is supported by a DOE OJI Award, a Sloan Foundation Fellowship, and by NSF 
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